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Abstract 

Both gravitational and mass induced neutrino oscillation mech- 
anisms provide possible resolutions to the Solar Neutrino Problem. 
The distinguishing feature between the two mechanisms is their de- 
pendence on the neutrino energy. We investigate the implications of 
this by computing the 8 B neutrino spectrum as determined from each 
mechanism using a realistic three-flavor evolution model. We find that 
in the limit of small #13 mixing angle, the differences are significant 
enough to observe in future solar neutrino experiments. 
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The Solar Neutrino Problem (SNP) has perplexed both astrophysicists 
and particle physicists for upwards of thirty years now. The measured flux 
of z/ e s incident on the Earth-bound detectors 0, [3|, [| remains in conflict 
with the standard solar model prediction (see || for a review). A plausible 
expanation of this discrepancy is that oscillations take place between electron 
neutrinos and neutrinos of differing species, thereby reducing the expected 
v e flux to an empirically acceptable value. 

At present there are two qualitatively distinct mechanisms which could 
give rise to such oscillations. One is the well-known Mikayev-Smirnov- 
Wolfenstein (MSW) mechanism [§, |7|, §|] which postulates that neutrinos pos- 
sess a non-trivial mass eigenbasis, in contrast to the assumptions of the Min- 
imal Standard Model. In this mechanism electron neutrinos produced at the 
core of the sun will, under certain conditions dependent upon the solar elec- 
tron density, undergo a resonance with other species of neutrinos whose flux 
then has no effect on earth-based detectors. Another mechanism proposed 
more recently f| hypothesizes that neutrinos posess a flavor-dependent cou- 
pling to the external gravitational field. This mechanism (recently dubbed 
the VEP mechanism [[TIJ) violates the Einstein Equivalence Principle (EEP), 
since it requires Gi — (1 + fi)G, where G is Newton's constant, i a flavor in- 
dex, and fi dimensionless parameters which characterize the degree of EEP 
violation, with each /j « 1. To ensure the full effect of three flavors, we 
must have fi ^ fj,i ^ j. For first generation neutrinos, we define fx = 0, 
i.e. G± = G. The VEP mechanism does not require neutrinos to have a 
non-degenerate mass-matrix. 

Both oscillation mechanisms rest on the assumption that the two neutrino 
eigenstates, flavor (\v)w) and mass/gravitational (|^)m,g), are related by an 
SU(Ng) transformation, 

W)w = H W)M,G , (1) 

where N g = 3 for three flavors. These states evolve according to the equation 
1 

ij- \ v)m,g = H m ,g W)m,g 
if r \v) w = {vjH M;G V 3 + A(r)} \u) w . (2) 

Flavor oscillations and resonances arise due to the off-diagonality of the mod- 
ified Hamiltonian V^H M)G V 3 + A(r), where A(r) = diag(V2G F N E (r),0,0) 
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is the term in the Hamiltonian corresponding to v e -e (i.e. charged-current) 
electroweak interactions. If we assume no CP-violation in the neutrino sec- 
tor, then the four-parameter matrix V3 reduces to a real, orthogonal rotation 
with three vacuum mixing angles 6> 12 , # 13 , and 6*23. 

The main difference between the MSW and VEP mechanisms manifests 
itself in the energy dependence of the evolution equations in (|2|). which 
respectively are 

H M = ^ I m\ I , (3) 

m\ 

H G = 2E\<j>\ I f 2 I , (4) 
h 

where a factor of unity has been subtracted out in each of (f|) and as it 
contributes only an overall unobservable phase. We have taken this factor 
to be f- (H Mt c)n, leaving the dynamics of the mechanism dependent of the 
eigenvalue differences Am^ = mf — m 2 , A/y = /, — fj. 

In this paper we consider the effect of these differing energy dependences 
on the suppression of the 8 B neutrino flux in the three-generation scheme. 
We show that for a small v e — > v T mixing angle #13 , there is a noticeable 
energy dependence in the shape and size of the suppression curves for MSW 
and VEP. This could be easily detected in present-day and future water 
detectors. For large #13, the suppression is energy-independent, and hence 
there is little variation in the spectrum. We shall take = 0©(r) in (|4|), 
where 0©(r) is the solar gravitational potential 1 . 

In order to study the suppression of z^ e s, we must determine the model- 
dependent survival probability for these neutrinos as they travel from the 
solar center to the Earth-based detectors. Averaged over 1 AU, the solution 



to (M) can be written HT2 



(p{ Ve - u e )) = £ \(v 3 ) u \ 2 {(PlzU 2 m m 



1 We do not consider the effects of ^-regeneration in the Earth (the "day-night effect"), 
nor do we consider the minimal contribution of the local supercluster. For discussions as 
to why these can be excluded, see pO[i or EH- 
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{(l-Pl)c? 2 C? 3 + P lS ? 2 4} 

+4i2<4i3 {A(l " + (1 - Pi)(l - ft)a?a4. + p 2sls} 

+4l3 {^^44 + ^(1 " A)4^3 + (1 - • 



The terms 



Sij = sin , Cjj = cos 9ij (6) 
refer to the vacuum mixing angles, while s m j,-, c m j 3 - are the analogous matter- 



enhanced trigonometric functions [14|. The referenced work expresses the 



functions in terms of MSW parameters; for VEP, we make the global substi- 
tution 

The functions Pi are the Landau-Zener jump probabilities [12| for non- 
adiabatic state transitions, and are inherently energy-dependent. 

The matter-enhanced mixing angles 9^ |14| are dependent upon the point 
of production of the solar u e . If this is above the resonance density (N e ) r ^ s 
for the specific 12- or 13-flavor transition (i.e. v e — > or v e — > u T ), then 
®u ~ y f • Hence, the probability in (H) reduces to 

(P(u e - V e )) = PlP2C 2 12 c\z + ^2(1 " Pl)s\Az + (1 " ^2)S?3 , (8) 

For large #13, the jump probabilities P2 vanish (i.e. adiabatic approximation 
for 13-transitions), and we are left with 

(P{"e - f.)> = 4 , (9) 

which is energy-independent (and, interestingly enough, independent of the 
12-transitions). Conversely, the small #13 limit of (Q) is 

(P(u e -> I/ e )) = (10) 

which shows distinct energy dependence on both 12- and 13-flavor transi- 
tions. By studying the two limits in question, it is possible to see exactly 
how the third flavor affects the dynamics of the oscillation mechanism. The 
small #i3 solutions will in some respects approximate the two-flavor mecha- 
nism (which is fully recovered in the limit #13 — > 0). At the other extreme, 



(5) 
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it is possible that the large #13 solutions can yield spectral distortions which 
match observed 8 B fluxes, but cannot be obtained via the two-flavor model. 

Since it was stated that the main difference between the VEP and MSW 
oscillation mechanisms should be visible in a study of the energy-dependent 
suppression of the neutrino flux, we show that striking differences can be 
observed for certain input parameters of the two models. 

Although the 8 B decay is one of the rarest nuclear reactions in the Sun, 
the resulting 8 B neutrinos are the easiest to study, since they have the widest 
spectrum of E G [0, 15]MeV, and are of sufficiently high energy to be observed 
by the v e -e scattering detectors (e.g. Kamiokande II, Superkamiokande, 
SNO). Hereafter, we provide a comparison of the 8 B flux curves as reduced by 
VEP and MSW, for both natural (A/31 > A/21, Am^ > Am^) and broken 
eigenvalue-hierarchies (A/ 31 < A/21, Am^ < Am^). 

While chemical detectors can ascertain the rate of solar neutrinos which 
reach the Earth, the water-based detectors provide extra pieces of informa- 
tion which can help pin down parameters. An accurate measurement of the 
incident 8 B flux can tell us such things as which direction the neutrinos came 
from, exact arrival times (assuming they originate from the Sun) and, most 
importantly, show the energy-dependence of the suppression mechanism at 
work. Due to the fact that the Cerenkov detectors have relatively high en- 
ergy thresholds, neutrinos whose energies are below these thresholds will not 
be visible. For example, KII has a threshold energy of E t h = 9 MeV [|1^ 



while SNO will have a lower one of E t h ~ 5 MeV. Thus, these experiments 
will be particularly useful if the parameter sets are such that suppression is 
visible in the high energy portion of the 8 B spectrum. 

The (chemical detector) counting rate R a for solar neutrinos from reaction- 
type a (e.g. 8 B, 7 Be decays, or the pp chain, hereafter 8 B, 7 Be, and pp 
neutrinos) is calculated via 

R a = f R& drr 2 C(r) dE <fP(E)<r(E)(P{v e -> v e ))(r,E) , (11) 

JO J E m i n 

where o~(E) is the detector-material absorption cross-section for neutrinos, 
4> a (E) the unreduced neutrino flux Hl5| , and £ Q (r) the fractional neutrino 
production rate at radius r 0. Note that the size of the flux curve does not 
reflect the resulting counting rate, so much as does its shape. The counting 
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rate in eq.flTTD is essentially calculated as 

iT = ]>», (12) 

i 

and so the same R can be obtained for small fluxes as well as for large ones, 
since the cross-section <7j = o~(E) increases with increasing energy. It is the 
location of the unsuppressed curve on the energy-axis which determines this. 
We have performed a numerical integration of Eq. ( O ) using the SSM data 
mentioned in the previous paragraph. Figures |l]-|5| show both MSW- and 
VEP-reduced 8 B fluxes for various allowed counting rates, as compared with 
the SSM (unreduced) flux of $^ A/ = 5.8 x 10 6 cm~ 2 s -1 . All mass-differences 
Am^ are expressed in units of eV 2 . 

The effects of the MSW/VEP mechanisms on solar neutrino depletion 
hinge on the existence of a matter-induced flavor-conversion resonance be- 
tween different neutrinos (e.g. v e — > v^). There are three possible resonances 
for three flavors: v e — > u^, v e — > u T , — > v T (and vice versa). For solar 
neutrinos, we are only concerned with flavor conversion from v e — > u x , where 
u x is either of the other two flavors.0 Resonances can occur if the v e s are 
created at an electron density (N e ) cr > (iV e )^ s , with i = 2, 3 for the other 
two flavors, where 

(Ne)ir = 7 ^v / 2|0(r)|A/i?cos2^ (13) 

For MSW, simply make the substition in eq. (0), and replace mass-eigenstate 
vacuum mixing angles with gravitational ones. 

Conversely, if (N e ) cr < (iV e )^| s , the v e will never undergo resonance, and 
will propagate as if in vacuum. We consider here only the case (N e )^ s > 
(N e ) max , i.e. the v e resonance density exceeds the maximal solar electron 
density. We refer to this behavior according to the following: single resonance 
((N e ) cr > (N e )™ s for % = 2 or 3 only), and double resonance ((iV e ) cr > [N e )™ s 
for both i — 2, 3). 

We note a marked difference between small and large angle solutions. 
Essentially, the large-angle solutions show very little variation in the fluxes. 
This is consistent with the form of @, which show the ^-attenuation to be 
energy independent. The least variation is visible in the double-resonance 

2 We do not consider arbitrary sterile neutrinos in this analysis. 
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case (fig. [I]), while we see more of a discrepancy in the models for a single- 
resonance (for both natural- and broken-hierarchies, figs [^, [|). This differ- 
ence is mostly in the low-energy neutrinos, though, and so would be difficult 
to detect. The curves for large #12, small #13 show similar energy indepen- 
dence, and thus are not presented here. 

In contrast, the small-angle solution shows very different behaviour. Fig- 
ure |] depicts the double-resonance reduced 8 B neutrino flux for the two 
models, with quite surprising dissimilarities. The two models show opposite 
energy-dependent reduction: MSW suppresses low energy neutrinos, while 
VEP suppresses high energy onesQ. Furthermore, the broken hierarchy case 
for the small angle region (fig. |5|) shows other intriguing behavior. As men- 
tioned earlier, the counting rate is determined by the product of the cross- 
section a(E) and the flux 4>{E). Here we see an example of where both large 
and small fluxes can represent the same counting rate. Whereas the natural 
hierarchy attenuated low energy neutrinos, here we see that the situation is 
reversed: the large curve shows the low-energy neutrinos largely uneffected 
by the MSW mechanism, while VEP leaves high-energy neutrinos alone. 

This radically different spectral distortion between the two models can in 
part be explained by the fact that (especially for small mixing angles) the 
adiabatic and non-adiabatic transitions are reversed. It is noted in [|l(| that 
the adiabaticity condition is violated for low energy neutrinos in VEP, while 
it is violated for high energy neutrinos in MSW. This is clearly reflected in 
the small #13 solutions of figs. f|. f| 

In summary, the inclusion of a third flavor in each mechanism yields a 
widely varying range of possible flux curves depending upon the values of 
the mixing angle parameters. The small angle solution offers variations in 
the structure of (P{v e — > Ve)), while still approximating the two-flavor limit 
(#13 — » 0). In the case of the small 13-mixing region for double flavor- 
resonances, the flux suppression is opposite for the two mechanisms. For 
an unbroken hierarchy, the VEP mechanism attenuates low energy neutrinos 
whereas MSW attenuates the higher energy range; in the case of a broken 
hierarchy the roles played by each mechanism are interchanged. The large 
#13 solutions are energy independent, thus making it difficult to distinguish 
between MSW and VEP in the high-energy portion of the spectrum, as most 

3 The rough nature of the curves is attributed to numerical variations, and most likely 
not a physical behavior. 
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spectral distortions occur here toward the lower end. Additional information 
from atmospheric observations [|16] and laboratory experiments [17] will then 
be essential in determining the underlying oscillation mechanism. 
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Figure Captions All captions list parameters used to obtain the re- 
duced counting rate R. Neutrino masses are expressed in units of eV 2 . 

Fig. 1: A/ 21 = 1(T 15 , A/31 = 4 x 1(T 14 ; Am 2 2l = 6.15 x 1(T 6 , Am 2 31 = 
5 x 1(T 5 ; 

Fig. 2: A/21 = , A/31 = 4 x KT 14 ; Am 2 21 = 6.15 x 1(T 6 , Am 2 31 = 

5 x 1(T 5 ; 

Fig. 3: A/ 21 = 1(T 8 , A/31 = 3 x 1(T 14 ;Am^ = KT 1 ,Am 2 31 = 

6 x 1(T 5 ; 

Fig. 4: A/21 = 2.7 x 1(T 15 , A/ 31 = 10" 14 ; Am 2 21 = 10" 5 , Am 2 31 = 
1.6 x 10" 4 ; 

Fig. 5: A/21 = 1<T 8 , A/31 = 1.31 x 1(T 14 ; Am 2 21 = 10" 1 , Am 2 31 = 
1.31 x 10~ 4 ; 
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(1e-15,4e-14),(8.15e-6,5e-5); R = 2.8 SNU, s(2th12) A 2 = 0.8, s13 A 2 = 0.4 
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Figure 2: 
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(1e-8,3e-14),(1e-1,8e-5); R = 1.60 SNU, s(2th12) A 2 = 0.8, s13 A 2 = 0.4 
800000 i 1 1 1 1 1 1 




160 



E (MeV) 



Figure 3: 
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(2.7e-15,1e-14),(1e-5,1.6e-4); R = 2.00 SNU, s(2th12) A 2 = 5e-3, s13 A 2 = 1e-3 
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Figure 4: 
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(2.7e-15,1e-14),(1e-5,1.6e-4); R = 2.00 SNU, s(2th12) A 2 = 5e-3, s13 A 2 = 1e-3 
800000 | 1 1 1 1 1 1 r- 
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Figure 5: 
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